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Abstract 

A tool for the identification of the shape of quantum dots is developed. By preparing a two- 
electron quantum dot, the response of the low-lying excited states to a homogeneous magnetic field, 
i.e. their spin and parity oscillations, is studied for a large variety of dot shapes. For any geometric 
configuration of the confinement we encounter characteristic spin singlet - triplet crossovers. The 
magnetization is shown to be a complementary tool for probing the shape of the dot. 

PACS numbers: 73.21.La,75.75.+a,71.70.Ej 
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I. INTRODUCTION 



The electronic structure of few-electron quantum dots (QDs) has been probed experi- 
mentally and provided us with new rich physics^&i. The basic physical picture consists 
of electrons effectively trapped by an in-plane parabolic artificial potential, a model which 
is compatible with the insensibility of the QDs to far-infrared radiation and is reflected in 
theory by the generalized Kohn theorem 5 * 6 -. Furthermore, for circular symmetric confine- 
ment a wide range of experimentally probed phenomena, such as addition energy spectr aM 
or inelastic light scattering'' measurements, could be explained by adopting such a model. 
The response of QDs to homogeneous magnetic fields has triggered rich scenarios. Among 
others, these are the magic angular momentum numbers, the configurations and the oscilla- 
tions of the total spin* 1 * as well as the parity oscillations for the ground state of a QD 1 ^ 11 
with changing field strength. The electron-electron interaction plays hereby a crucial role. 

Introduction of the anisotropy breaks the rotational symmetry thereby altering the dot's 
electronic and dynamical properties. The shell structure in the addition energy spectr a 2ll2lls 
and the level clustering in the two-electron QD excitation spectru m 14 ! 15 are eliminated and 
dynamically the anisotropy serves as a rapid path to chao a 14 i 15 . The spin configurations 
obey the Hund's rule for small QDs while for larger ones they change to a more Pauli- 
like behavio r 16 ! 17 . In three dimensions the anisotropy has been introduced with a magnetic 
field which therefore controls the symmetries^*. However, despite the enviable advances in 
revealing* and explaining 2 - the electronic structure of parabolic QDs, the question for the 
exact shape of a dot, as it is realized experimentally, can not be answered precisely. A 
major approximation often made, suggests an analogy between the geometric configuration 
of the device and the effective confinement potential 2 *!*. However, it is being admitted that 
it is nearly impossible to derive the exact theoretical model of the confinement from the 
geometry of the setup of the device. Even within circular geometry the precise electronic 
properties depend on the details of the device*. 

In the present article the properties of two-electron parabolic QDs in magnetic fields are 
investigated for the full deformation regime from circularly shaped to wirelike dots. In Sec. 
II we provide the Hamiltonian within the framework of the effective mass approximation and 
describe its general symmetries as well as our computational approach. In Sec. Ill we present 
and discuss our results. In particular, the properties of the two-electron anisotropic QD (Sec. 
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III. A) and the isotropic two-electron QD in magnetic field (Sec. III.B) are summarized in 
order to become acquainted with the indivudual effects of the anisotropy and magnetic 
field in our two-electron system. Sec. III.C provides the response of the spin and parity 
symmetries of the ground and few low-lying states as well as the magnetization to the 
magnetic field for the full deformation regime addressed here. We encounter unique features 
for every dot shape. Ultimately, we develop a tool for the diagnosis of the shape of dots, as 
they are realized experimentally. 

II. HAMILTON! AN AND COMPUTATIONAL APPROACH 

In the effective mass approximation the conduction band electrons confined in a two- 
dimensional general parabolic quantum dot in a magnetic field B = (0, 0, B) are described 
by the Hamiltonian H = Hcm + 7~tr with 

Hcm = -— (P + 2eA(R)) 2 + m e u 2 ( cos 2 : X 2 + sin 2 F 2 ) (1) 
4m e 

H r = ^-( P + ^A(r)) 2 + ^c 2 (co S 2 x 2 + S m 2 y 2 ) + (2) 

where we choose for the vector potential the symmetric gauge A(r) = | (B x r) and 
e,m e ,e,u ,(p are the electron charge, effective mass, dielectric constant, the characteristic 
frequency and the anisotropy parameter, respectively. Small and capital letters refer to 
the relative and center of mass degrees of freedom, respectively. In the following we adopt 
the typical values for a GaAs dot (effective Bohr radius a* B = 9.8 nm, effective Hartree 
Ha* = 11.8 meV, hu Q = 4.96 meV and 1 effective unit (e.u.) of field strength corresponds to 
6.925 Tesla). Quantization of Hcm is straightforward. Direct observation of the electronic 
properties due to H r via far infrared spectroscopy is prohibited, since radiation in the dipole 
approximation contains only CM degrees of freedom. In the following we focus on the 
non-trivial relative motion H r - Parity (r — > — r) and spin are interrelated symmetries due 
to the Pauli exclusion principle and we encounter spin singlet eigenf unctions with even 
spatial symmetry \l/(r) = \&(— r) and spin triplet eigenf unctions with odd spatial symmetry 
^(r) = — r). The two characteristic frequencies of the confinement are u x = uo Q coscf) 
and ujy = uj sincf). For = 45° the dot has a circular shape. With increasing it obtains 
an elliptic shape and ends in a wirelike dot for — > 90° {uj x — > 0, u y — > uj ). 
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To investigate the two-electron QD we solve the corresponding Schrodinger equation using 
a full configuration interaction approach with the harmonic oscillator-like basis set 

®n x n y = A{n x ,n y )H nx {JF x x)H nv {^y) x ^ « 2+ « x ~% >» (3) 

which leads to an algebraic eigenvalue problem. In Eq. (j3J) A(n x ,n y ) is the normaliza- 
tion constant, c\ = M\U)\/c, C3 = M2UJ2/C, c 2 = 2iiM\oj\M 2 ijj<ii 'c, c = /i 2 MxM^ojxoj^ + 1, 
H = -2L/(m e u p), A = [m e uj L{2 + L 2 )]/[4(cos(20) - p)\ , M h2 = m e p/{p - cos{2<p) =p L 2 ), 
^1,2 = (c^o / \/2) yl + L 2 =(= p, L = eB/m e ujQ and p = + £ 2 ) 2 — sin 2 {2(f)). All the units 
are scaled appropriately. The eigenf unctions \t„ litl2 of the charged anisotropic harmonic 
oscillator (corresponding to 7i r without Coulomb repulsion) can be obtained analytically^ 2 -. 
The exponential function of our basis set Q represents exactly the exponential part con- 
tained in ty ni>n2 . The Hermite polynomials contained in ^ ni ,n 2 can be obtained through a 
finite superposition of the corresponding Hermite polynomials in the orbitals <& nxjny . In or- 
der to arrive at the eigenstate ^ ni ,n 2 we need to superimpose Hermite polynomials of degree 
n x = n y = rii + n 2 for both x and y directions. For the evaluation of the Hamiltonian matrix 
elements (without the Coulomb interaction) we subtracted from Eq. (J2J the Hamiltonian 
whose eigenfunctions are The remaining electron-electron integral was evaluated 

numerically. The Coulomb repulsion term was replaced by an auxiliary Gaussian integral 
and the resulting three-dimensional integral was evaluated by combining Gauss-Hermite ex- 
act quadratures with respect to x and y directions and a Gauss-Kronrod quadrature in the 
remaining auxiliary variable. For details we refer the reader to Refs. 1 1 4JJ 1 5l . 



III. RESULTS AND DISCUSSION 
A. No magnetic field 

Before we proceed with the analysis of our results let us briefly discuss some principal 
properties of our system without magnetic field illustrated in Fig. 1 (see also Refs. 

□a. 

For cf) = 45°, Ti. r is rotationally symmetric and L z is conserved. The sequence of excited states 
with increasing energy is: (m; S) = (0; 0),(±1; 1), (±2; 0),(0; 0),(±3; 1),(±1; 1),. . ., where m 
and S are the magnetic quantum number and the total spin, respectively. Introduction of 
the anisotropy splits the doubly degenerate levels (for m 7^ 0) and leads to spin singlet - 
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FIG. 1: Schematic diagram of the energy of the ground and first seven excited states for B = 
with increasing (j). The vertical line indicates the angle for which the system is integrable. 

spin triplet (ST) crossings. At uj v : uj x = 2 : 1 the system becomes integrable due to the 
non-linear constant of motion 



A = {L g , Px } + t fyx 2 - J— (4) 

The operator A commutes with the x-parity and anticommutes with the ^/-parity. These 
symmetries result in ST degeneracies which occur systematically in the excitation spectrum 
(see Fig. 1). For — > 90° the low-lying states arrange themselves in energetically well- 
separated pairs involving one spin singlet and one spin triplet state. 



B. Isotropic case in a magnetic field 

The two-electron isotropic QD in a magnetic field has been well studie d 1 ^ 1 ?? 2 ?! 21 and 
particular analytical solutions have been found^. The system is integrable, with L z being a 
constant of motion. A major feature of this system are the parity oscillations of the ground 
state 1 ' , being a clear imprint of the Coulomb interaction. By increasing the magnetic field, 
confinement strengthens thereby increasing the Coulomb energy between the two electrons. 
Consequently, the system attempts to establish itself in a lower energy by pushing apart the 
electrons into states with higher \m\. Hence, by increasing the field the ground state pumps 
itself into an energetically favourable state, thereby flipping the spin, with the following 
sequence: (m; S) = (0; 0),(— 1; 1),(— 2; 0),(— 3; 1),. . ., i.e. we encounter ST oscillations for 
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FIG. 2: (Color online) Domains of spin multiplicity in the (B,4>) plane for (a) the ground state 
and (b) the fourth excited state. Color indicates the energy difference AEq and AE4, respectively, 
on a logarithmic scale. Deep red (dark regions) and yellow to green (bright regions) correspond 
to large and small spacings, respectively. The yellow (bright) curves form the borders between the 
different ST-symmetry domains. 

the ground state. For excited states an increasing number of crossings occur— due to the 
participation of states covering an increasing range of m. For increasing field strength, the 
level spacing AEi = E i+1 — Ei (i determines the degree of excitation and takes even values 
% = 0,2,4, ... within our study) oscillates between zero and a maximum value AE max . For 
stronger fields the amplitude AE max reduces and the slope of the oscillating energy curves 
converges to a constant value. This fact leaves its fingerprints in the magnetization which 
for zero temperature is defined by M(B) = — (7^) with E (B) being the ground state 
energy. The ST crossovers will then appear as steps in the magnetization. With increasing 
temperature these steps smear out and result in a smooth diamagnetic behavior. Typically, 
for T < 0.1K the magnetization's signal reveals the ST oscillations. 

C. Anisotropic configuration in the magnetic field 

Let us now proceed with the introduction of the anisotropy in the two-electron isotropic 
QD in the magnetic field. The rotational symmetry breaks and a large number of avoided 
crossings between states with identical symmetry occur. Fig. 2(a) shows the spin multiplicity 
S = 0, 1 of the ground state in the (B, <fi) plane. We observe alternating islands of different 



6 



total spin. For = 45° and the given range of magnetic field we receive ST oscillations 



in correspondence with Ref. 



1CM . Introducing the anisotropy, for B = 0, splits the first 



excited doubly degenerate triplet state (£1,2) ■ As a result AE decreases gradually with 
and in the wirelike limit — > 90°, E\ approaches E (A-Eq ~~ * 0) forming an energetically 
well-separated ST pair (see Fig. 1). For fields B < 1 the border curves separating the 
domains of different spin multiplicity depend very weakly on the field strength. For higher 
field strengths 1 < B < 2 a more significant change can be observed as e.g. the third (5 = 0) 
and fourth (5 = 1) domain slightly spread and higher domains such as the fifth (5 = 0) 
domain for > 54° (see upper left corner of Fig. 2(a)), are suppressed. 

The third and fourth excited state for = 45° and B = form a doubly degenerate spin 
singlet pair (-£3,4). By introducing the anisotropy the corresponding energy splits and at 
uj y : uj x = 2 : 1 the energy E% becomes degenerate with the energy E 2 of the lower excited 
spin triplet state, an effect of the constant of motion and its associated symmetries for this 
particular anisotropic configuration (see Sec. III. A). Hence, for the second excited state 
(corresponding to E 2 ), we observe a ST transition already in the absence of the field, being 
a pure consequence of the geometric change of the confinement. For stronger anisotropies 
(^ > 2) the energies E 2 , E 3 of the two states converge again and form an energetically well- 
separated pair of energies. By introducing the magnetic field, the spin-symmetry domains 
change more significantly compared to the ground state. Hence, for E 2 , the first (5 = 1) 
domain, with increasing 0, will be suppressed at u y : u x = 2 : 1 by the following (5 = 0) 
domain. The domains, for stronger fields, slightly widen after the suppression of the first 
5=1 domain. 

Fig. 2(b) shows the domains of spin multiplicity 5 — 0, 1 in the (B, 0) plane for the 
fourth excited state. The fourth excited state (£4) is a spin singlet, for B = and = 45°, 
and at u y : u x = 2 : 1 (0 w 63.4°) it crosses with the energy curve E 5 (see Figs. l,2(b)) 
thereby becoming a spin triplet. For stronger anisotropies and due to the large number of 
excited states, an 'accidental' crossing with a spin singlet state occurs and the eigenstate 
belonging to E\ restores its initial parity and spin-multiplicity. In the wirelike limit, a third 
pair of energetically well-separated states is formed. For weak magnetic fields and changing 
0, the first (S = 0) domain disappears at u y : u x « 2 : 1 and a (S = 1) domain follows. For 
> 75° the 5 = 1 domain is suppressed, due to the above-mentioned 'accidental' crossing, 
by a second 5 = domain. For stronger fields, the domains evolve smoothly and slightly 
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FIG. 3: The magnetization M(B) for various anisotropics (ft. The different curves are offset by 0.3. 

widen after the accidental crossing at ~ 75°. From Fig. 2(b) it can be seen that the border 
lines separating the spin domains in the (B, (ft) plane strongly depend on the anisotropy (ft. 

We therefore conclude that the map of ST domains in the (B, (ft) plane for low-lying states 
shows unique features that are characteristic for the shape of the dot. 

An important tool to study the response of the QD to the presence of the external 
magnetic field is its magnetization M. Calculations within the Hartree approximation for 
three and four electrons have shown that M is sensitive to the number of electrons and the 
shape of the dot>2i Fig. 3 shows M(B) for the two-electron QD and various anisotropics. 
For (ft = 45° the curve shows four steps whose positions correspond to the ground state ST 
crossovers. For lower magnetic fields the steps are more pronounced due to the very different 
slopes of the crossing energy curves. Hence, the magnetization contains also information on 
AE (B). By introducing the anisotropy, AE (B = 0, (ft) decreases and the amplitude of the 
oscillations AE max also decreases thereby reducing the height of the steps. Consequently, the 
steps with increasing (ft gradually vanish and finally, we obtain a smooth curve for (ft = 85°. 
The magnetization, therefore, can reveal also significant information on the shape of the dot, 
requiring the study of the ground state only. 

The above discussion holds if the contribution of the Zeeman term Eg(B) = g*fisBS z , 
(with hb being the Bohr magneton and g* = —0.44 for GaAs) is neglected. E$ splits the 
threefold degeneracy of the spin triplet states while the spin singlet states remain unchanged. 
Hence, we expect a suppression of the spin singlet ground states for field strengths where 
AEq(B) is smaller than Es(B). As a result, for high magnetic fields, E$ dominates and the 
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FIG. 4: (Color online) Domains of spin multiplicity of the ground state in the (-B,0) plane at the 
presence of the Zeeman spin splitting. 

ground state is always a spin triplet^. Fig. 4 shows the spin multiplicity of the ground state 
in the (-8,0) plane in the presence of the Zeeman term. For = 45° the ground state shows 
ST oscillations. Starting from B = and increasing B, the interval of field strengths of the 
first S = domain is approximately the same as for g* = 0. However, the following S = 1 
domains are enlarged and the second S — domain shrinks significantly. Introduction of 
the anisotropy causes a reduction of AE which is not noticeable for low field strengths but 
leads, for stronger fields, to the disappearance of the second spin singlet island already for 
~ 48°. Hence, apart from this spin singlet island, the ground state is given by a spin 
triplet for B > 0.3, in the regime from circular to intermediate anisotropies (0 < 65°). For 
higher anisotropies, AE decreases significantly, resulting in a smooth elimination of the 
first spin singlet domain as approaches the wirelike limit. 

An important issue to be discussed is the effect of the Zeeman splitting for higher excited 
states in the presence of anisotropy. Clearly, for larger fields, the ST oscillations will be sup- 
pressed in favour of the spin triplet symmetry due to the dominance of the Zeeman splitting. 
For lower fields, however, the ST oscillations persist with varying 0. This statement holds 
typically for B < 1 and a significant number of excitations. The results obtained above for 
g* = therefore transfer immediately to the case g* = —0.44. 



IV. CONCLUSIONS 



To conclude, a tool for the identification of the shape of a QD has been developed, by 
preparing a two-electron QD and studying its response, when switching on a homogeneous 
magnetic field. We examined low-lying excited states, that are accessible experimentally, 
with e.g. tunneling spectroscopy. By applying a bias voltage between the two reservoirs 
weakly coupled to the dot, the Coulomb blockade region can become narrower and the extra 
tunneling channels, corresponding to the excited states, can be detected as an increase in 
the current 4 . The parity and spin symmetries of these states change uniquely with varying 
the shape 0, for field strengths up to several Tesla. Hence, the detection of the positions 
of the ST crossovers for low-lying excited states allows us to conclude uniquely on the 
shape of the QD. As an alternative, independent possibility to confirm the results we 
suggest magnetization measurements. Our investigation yields the 0-dependence of the 
magnetization showing a gradual transition from a step-like to smooth behavior. In the 
presence of the Zeeman term (g* = —0.44) and for weak to intermediate field strengths our 
results remain practically unchanged and can be equally applied. In the strongly anisotropic 
regime, already the ground state ST crossover curve in the (-8,0) plane shows a monotonic 
change with 0. The method developed here applies, therefore, to the full deformation regime. 
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